N PS  ARCHIVE 
1962 
SALLADA,  W. 


RAPID  APPROXIMATION  OF 
SERVOMECHANISM  FREQUENCY  RESPONSE 

WILLIAM   F.  SALLADA 


VAUGHN  E.  WILSONf  JR. 


LIBRARY 

m  tm  fv  kmhv  i  lao a»v  0  *•  naval  postgraduate  school 

UUOLtY  KNOX  LIBRARY  wo-  ,RI  nA 

NAVAL  POSTGRADUATE  SCHO 
MONTEREY  CA   93943-5101 


RAPID  APPROXIMATION  OF  SERVOMECHANISM 
FREQUENCY  RESPONSE 

William  F.  Sallada 
and 
Vaughn  E.  Wilson,  Jr. 


RAPID  APPROXIMATION  OF  SERVOMECHANISM 
FREQUENCY  RESPONSE 


by 

William  F.  Sallada 
Lieutenant  Commander,  United  States  Navy 

and 

Vaughn  E.  Wilson,  Jr. 
Lieutenant,  United  States  Navy 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 
ELECTRICAL  ENGINEERING 


United  States  Naval  Postgraduate  School 
Monterey,  California 

1962 


LIBRARY 
U.S.  NAVAL  POSTGRADUATE  SCHOOL 
'  MONTEREY,  CALIFORNIA 
RAPID  APPROXIMATION  OF  SERVOMECHANISM 

FREQUENCY  RESPONSE 
by 

William  F.  Sallada 
and 
Vaughn  E.  Wilson,  Jr. 

This  work  is  accepted  as  fulfilling 
the  thesis  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 
ELECTRICAL  ENGINEERING 

from  the 
United  States  Naval  Postgraduate  School 


ABSTRACT 

The  purpose  of  this  thesis  was  to  develop  a  technique  by  which  the 
frequency  response  of  a  servomechanism  might  be  rapidly  evaluated  by 
graphical  methods.   The  approach  consisted  of  developing  several  frequency 
response  parameters  related  to  the  location  of  the  poles  and  zeros  of  the 
closed  loop  transfer  function  of  a  servomechanism.   The  systems  considered 
were  the  pure  second  order  (one  pair  of  complex  conjugate  poles),  second 
order  with  one  zero  (one  pair  of  complex  conjugate  poles  and  one  real 
zero) ,  and  a  pure  third  order  (one  pair  of  complex  conjugate  poles  and 
one  real  pole).   For  the  latter  two  systems,  the  complex  poles  and  the 
frequency  were  redefined  as  referred  values  relative  to  the  real  pole 
or  zero. 

Loci  of  constant  frequency  response  parameters  were  mapped  onto  a 
complex  plane;  or  loci  of  constant  imaginary  referred  pole  locations  were 
mapped  onto  an  auxiliary  complex  plane,  depending  upon  the  specific 
system  considered.   Using  the  points  obtained  from  the  appropriate 
relationship,  the  frequency  response,  M  and  N  contours,  could  be  rapidly 
and  easily  sketched.   The  M  and  N  contours  could  then  be  used  as  a  basis 
for  rapid  system  evaluation  in  the  frequency  domain. 

The  authors  wish  to  thank  Professor  George  J.  Thaler  for  bringing 
to  their  attention  the  thesis  by  Major  Walter  B.  Patton,  USMC ,  and 
LT  William  B.  Abbott,  III,  USN,  without  which  this  thesis  topic  would 
probably  not  have  been  considered  /l/. 

The  authors  are  deeply  indebted  to  Assistant  Professor  Robert  D.  Strum, 
who  acted  as  their  thesis  advisor,  for  the  friendly  support,  encour- 
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1.   INTRODUCTION 

Each  year  the  importance  and  utilization  of  servomechanisms  and 
automatic  control  systems  increase.   It  naturally  follows  that  more 
papers  are  submitted  pertaining  to  the  design,  analysis,  synthesis, 
and  compensation  of  servomechanisms  and  automatic  control  systems. 
Many  of  these  papers  and  existing  texts  discuss  details  and  laborious 
calculations  which  are  certainly  a  prerequisite  to  analyze  or  synthesize 
the  operation  of  a  servomechanism.   However,  in  many  cases,  a  rapid 
approximation  of  system  operation  would  be  beneficial,  and  possibly 
only  an  approximation  would  be  required  upon  which  an  intelligent 
estimate  could  be  made  on  the  performance  of  the  system  under  study. 
There  are  only  a  few  papers  known  to  the  authors  in  which  the  use  of 
approximations  are  presented  from  which  an  estimate  of  system  performance 
can  be  made,  e.g.,  Floyd's  method  111.     One  such  paper  by  Patton  and 
Abbott  III   presents  a  method  for  a  rapid  approximation  of  servomechanism 
transient  response. 

A  corollary  should  exist  to  the  Pat ton- Abbott  development  in  which 
a  rapid  approximation  of  the  frequency  response  of  a  servo  could  be 
determined.  The  use  of  frequency  response  analysis  is  still  an  important 
criterion  in  measuring  servo  performance,  although  the  transient  response 
approach  has  become  very  popular  and  widely  used  in  recent  years,  almost 
to  the  point  of  neglecting  frequency  response.  There  are  many  ways  of 
determining  frequency  response  from  transient  response,  but  lengthy 
calculations  are  necessary. 

It  Is  the  purpose  of  this  paper  to  present  a  graphical  means  for 
determining  a  rapid  approximation  of  servomechanism  frequency  response. 
This  approach  eliminates  the  use  of  the  Bode  diagram  and  Nichols' 
chart.  The  M  and  n  contours  will  be  obtained  directly  from  the  closed 
loop  transfer  function  of  the  servo. 


This  development  will  parallel  the  Patton-Abbott  approach  so  that 
from  the  same  given  data  (closed  loop  transfer  function,  factored  form) 
it  is  possible  to  obtain  rapidly  the  transient  response  from  their 
paper  and  the  frequency  response  from  this  paper.   This  paper  will 
discuss  in  sufficient  detail  the  development  of  the  frequency  response 
approach  so  that  no  reference  to  the  Patton-Abbott  paper  will  be  required, 
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2    GENERAL  DEVELOPMENT 

Since  the  time  domain  or  transient  response  is  used  so  widely  in 
the  initial  phases  of  servo  study,  all  developments  in  this  paper 
originated  with  the  closed  loop  transfer  function  expressed  in  Laplace 
notation  (s-plane)  with  all  initial  conditions  equal  to  zero.   The 
transformed  equation  was  then  converted  into  frequency  response  notation. 
For  this  to  be  a  valid  conversion,  certain  basic  definitions  are  in  order, 
plus  some  additional  restrictions  and  assumptions  used  in  this  devel- 
opment.  The  frequency  response  of  a  system  was  defined  in  part  by 
considering  only  steady  state  conditions.   Thus,  all  transient  conditions 
have  vanished.   To  convert  the  transfer  function  from  Laplace  notation 
into  the  frequency  response  equation,  it  was  only  necessary  to  make  a 
substitution  of  the  variable.   This  was  based  on  the  assumption  that 
the  input  was  a  pure  sinusoid  and  that  the  output  was  also  a  pure 
sinusoid  since  the  steady  state  condition  exists.   Furthermore,  this 
forces  the  restriction  that  the  system  must  be  linear.   The  complex 
variable  s  ■  o  +  jou  then  becomes  the  imaginary  frequency  variable 
since  the  real  part  equals  zero  (a  ■  0).   Therefore  s  ■  ju);  and  when 
this  variable  substitution  was  made  in  the  closed  loop  transfer 
function,  the  frequency  response  equation  resulted. 

2. 1  Basic  Definitions  and  Systems  Characteristics 

At  this  point,  a  tabulation  of  the  various  definitions,  restric- 
tions, and  assumptions  would  be  in  order: 

a.  The  system  was  linear. 

b.  The  system  was  stable. 

c.  The  system  was  characterized  by  unity  feedback,  i.e., 
M  ■  1  at  ou  ■  0. 

d„   The  closed  loop  transfer  function  was  known  in  its  factored 
form. 

Based  on  the  above,  some  general  frequency  response  characteristics 
could  be  stated. 


2.2  M-Contour 

a.   At  zero  frequency,  M  is  unity. 

b»   At  zero  frequency,  the  slope  of  the  M  contour  is  zero. 

c.  As  the  frequency  approaches  infinity,  M  approaches  zero. 

d.  As  the  frequency  approaches  infinity,  the  slope  of  the  M 
contour  approaches  zero. 

e.  At  the  maximum  or  minimum  value  of  M,  the  slope  of  the  M 
contour  is  zero. 

2.3  N-Contour 

a.  At  zero  frequency,  N  is  zero. 

b.  As  the  frequency  approaches  infinity,  N  approaches   "0       tt 
degrees. 

(m  is  the  number  of  zeros  and  n  is  the  number  of  poles  of  the  closed 
loop  transfer  function) 

2.4  Parameters 

To  display  graphically  the  frequency  response,  parameters  were 
selected  to  define  specific  points  on  the  M  and  N  contours.   The 
selected  parameters  were  defined  as: 

a.  M  :   The  value  of  M  where  the  slope  of  the  M  contour  was  zero 

P 

(Does  not  include  the  value  of  M  at  zero  frequency  when  M  ■  1.0). 

b.  ou  :   The  value  of  omega  (U)  ■  2nf)  where  the  slope  of  the  M 

P 

contour  was  zero.   (Does  not  include  to  ■  0). 

c.  N  :   The  value  of  N  at  (!)  . 

P  P 


d. 


t)W 


e.  cu. 

bw 

f.  N. 

bw 


The  bandwidth  point,  or  when  M  ■  .707. 

The  bandwidth  frequency  (when  M  ■  .707). 

The  value  of  N  at  the  bandwidth  frequency, 
g.  M.:   The  notation  when  M  ■  1.0,  except  at  (i)  ■  0. 
h.   U)  :.   The  value  of  u>  when  M  -  1.0,  except  at  go  «  0. 
i.   N-:   The  value  of  N  at  it),. 
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j.   id  :    The  undamped  natural  angular  frequency. 

k.   M   r,   The  value  of  M  at  the  undamped  natural  angular  frequency  (ou  ) 
(Dn  n 

1.   N   :   The  value  of  N  at  the  undamped  natural  angular  frequency  (id  ) 
uun  n 

The  above  selected  parameters  defined  the  general  shape  of  the  M 

and  N  contours.   However,  it  was  discovered  that  with  only  these  points, 

the  M  and  N  contours  were  in  some  cases  too  general;  the  contours  had 

insufficient  shape  to  be  of  significant  value.   Better  approximation  of 

the  M  and  N  contours  were  provided  by  the  selection  of  additional  M,  N, 

and  id  points  to  "fill-in"  the  shape  of  the  contours.   It  was  discovered 

that,  mathematically,  the  choice  of  cd  as  the  basis  for  these  additional 

P 
points  reduced  the  M  and  N  equations  to  relatively  simple  relationships 

and  insured  that  the  additional  points  were  located  in  the  critical 

regions  of  the  contours.   The  benefits  gained  by  the  additional  points 

more  than  compensate  for  the  slight  additional  labor  by  providing  a  much 

more  accurate  approximation  to  the  M  and  N  contours.   The  additional 

parameters  selected  were: 

m.   (D  :   An  arbitrary  frequency  proportional  to  CD   (CD  K  k  (l)  ), 

where  k  was  a  selected  constant  as  discussed  later. 


n,   M  :    The  value  of  M  at  CD  . 

0.  N,:    The  value  of  N  at  CD . 

A  sketch  of  the  general  M  and  N  contours  is  shown  in  Fig.  2-1. 

2 . 5   Classes  of  Servomechanisms  Studied 

From  Fig.  2-1,  it  can  be  seen  that  it  was  relatively  easy  to 
construct  an  approximate  frequency  response  if  the  above  defined  para- 
meters were  known   This  paper  contains  the  above  parameters  presented  in 
graphical  form.   This  paper  also  contains  the  technique  utilized  in  the 
parameter  determination  such  that  the  frequency  response  can  be  rapidly 
approximated  for  specific  classes  of  systems.   The  three  systems  covered 
in  this  paper  were: 

1.  Pure  Second  Order  -  The  closed  loop  transfer  function  consisted 
of  two  complex  poles  only. 


M 


U) 


Fig.  2-1  General  M  and  N  Contours  Showing  Parameters  Chosen 
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II.  Second  Order  with  One  Zero  -  The  closed  loop  transfer  function 
consisted  of  one  pair  of  complex  poles  and  one  real  zero. 

III.  Pure  Third  Order  -  The  closed  loop  transfer  function  consisted 
of  one  pair  of  complex  poles  and  one  real  pole. 

2  6   Presentation  of  Parameters 

The  graphical  presentations  for  the  pure  second  order  system  were, 
for  the  most  part,  loci  of  constant  values  of  the  various  parameters 
plotted  on  the  s- plane.   For  the  other  two  systems  covered  in  this 
paper,  it  was  necessary  to  express  the  coordinates  of  the  complex  poles 
as  a  function  of  the  magnitude  of  the  real  pole  or  zero.   In  accomplishing 
this,  the  angular  frequency,  (i),  also  became  a  function  of  the  real  pole 
or  zero  location-   Thus  the  magnitude  and  phase  shift  for  the  frequency 
response  was  shown  to  be  a  function  of  the  complex  pole  location  relative 
to  the  real  pole  or  zero.   The  graphical  presentations  for  the  second 
order  plus  one  zero  and  for  the  pure  third  order  systems  were  plotted  in 
two  ways.   First,  by  loci  of  constant  values  of  the  various  parameters 
plotted  on  a  referred  complex  plane,  that  is,  with  ordinate  and  abscissa 
being  the  complex  pole  location  relative  to  the  real  pole  or  zero. 
Second,  by  loci  of  constant  referred  imaginary  part  of  the  complex  pole 
plotted  against  referred  real  part  of  the  complex  pole  as  abscissa  versus 
the  chosen  parameter  as  ordinate. 

2 . 7   Digital  Computer  Utilization 

The  solutions  of  the  equations  that  will  be  presented  later  when 
a  specific  system  is  discussed  were  obtained  from  a  Control  Data  Corp- 
oration high-speed  digital  computer  Model  1604,  located  at  the  U.  S.  Naval 
Postgraduate  School,  Monterey,  California  (Appendix  A).   Without  the 
availability  of  a  high  speed  digital  computer,  it  is  doubtful  that 
this  paper  could  have  been  developed  without  the  utilization  of  an 
analog  computer  plus  excessive  labor. 
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3    THE  PURE  SECOND  ORDER  SYSTEM 

The  pure  second  order  system  treated  in  this  section  was  pre- 
viously defined  as  a  system  whose  closed  loop  transfer  function  consisted 
of  only  two  complex  conjugate  poles.  This  system  has  been  extensively 
treated  in  most  servo  textbooks,  usually  in  terms  of  the  damping  ratio 
and  the  undamped  natural  frequency..  Here,  however,  it  will  be  treated 
by  the  rectangular  location  of  the  closed  loop  poles.  The  closed  loop 
transfer  function  of  the  pure  second  order  system  was  stated  as  /!/: 

r         '      _____ 
■JS)  '  (v  '  (3-1) 

The  poles  of  this  system  were  plotted  on  the  s-plane  in  Fig.  3-1. 
Converting  Eq .  (3.1)  to  the  frequency  response  form,  it  became 


'  (  '  '    /-"       ,  J  *  '  (3.2) 


3. 1  M  and  uu  Calculations 

Equation  (3.2)  rearranged  became 


F(;u}=       *  +  k' 


,_ 


C3  3^ 


from  which 


/    = 


!  j      i 


(3.4) 


Note  that  if  the  system  gain  were  constant,  then  "a"  and  "b"  would  be 
fixed  and  the  magnitude  would  be  a.  function  only  of  system  frequency, 
i.e.,   M  =  f ( _) 


s -plane 


X-- 


a. 


|U-> 


y<y 


Fig.  3-1   Poles  of  the  Pure  Second  Order  System  (s-plane) 
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3.2  M  and  to  Determination 
_JE E 

to  was  obtained  by  finding  the  value  of  CO  where  the  slope  of  the 
P 
M  versus  to  curve  was  equal  to  zero.   For  mathematical  convenience,  this 

was  done  by  first  squaring  Eq .  (3.4),  then  taking  the  first  derivative 

and  setting  the  numerator  equal  to  zero. 


-  (a  +  b*")  [t  (a-h  bX-(J-)('X co)  +  8  clco  ]  —  0  (3.5) 

Simplifying,  this  becomes 

W  =  o>p  =-Vb*'-^a                       (3.6a) 

and   to  =  0  (see  Section  2.4b)                (3.6b) 

Substituting  Eq .  (3.6a)  in  Eq .  (3.4)  gave  the  value  of  M  at  (l) 

.  .      _  CL  +  & 


Note  that  from  Eq.  (3.6),  M  would  not  occur  whenever  "a"  was  less  than 
"b" ,  i.e.,  when  the  real  part  of  the  pole  was  greater  than  the  imaginary 
part.   (Figs.  3-2  and  3-3) 


3.3  cu.   Determination 

bw 

The  equation  for  the  bandwidth  frequency  was  obtained  by  substituting 
the  bandwidth  magnitude  (M  =   707)  into  Eq.  (3.4)  and  solving  for  the 
frequency   This  resulted  in 


iJOi 
fcw 


UAl4-\la4+b4  (3-8) 


(Fig.    3-4) 
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3.4      M        and  ou     Determination 
uun _n 

The   undamped   natural    frequency   was    defined   as 


ti->n~    \0         '  (3.9) 


This  equation  expressed  uu  directly  in  terms  of  "a"  and  "b" .   (Fig.  3-5) 
The  substitution  of  uo  into  Eq .  (3.4)  yielded 

■    •' 


M. 


(3.10) 


Equation    ( 3, 10)expressed  M        in    terms   of    the    closed    loop   poles   of    the 

uun 

system.   (Fig.  3-6) 


3  5   U)  Determination 

As  defined  previously,  u>  was  the  frequency,  other  than  zero,  at 
which  M  was  unity.  Letting  M  =  1.0  in  Eq .  (3.4)  and  solving  for  u> 
yielded 


(3.11) 


Simplifying  this  expression  by  recalling   OJp  ==T/jb*"—  dL~ 

resulted  in  Oj         '!        (  (3.12) 

(Fig-  3-7) 

3  6  M,  and  U),  Determination 
k      k 

As  stated  in  Section  2,  intermediate  values  of  M  and  N,  in  addition 

to  those  determined  above,  would  result  in  a  smoother  plot  approximating 

the  M  and  N  contours.   The  values  of  k  chosen  for  this  class  of  servo 

were  .5  and  1.5    These  particular  values  were  selected  as  being  those 
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points  which  would  be  most  beneficial  in  smoothing  the  plot  of  the 
M  and  N  contours^  and  because  the  VL     equation  was  the  same  for  k  =  .5 
as  tor  k  =  1.5.    In  Section  2, 


and  Uk~i%"p^iLSwf 

(Fig      3-7) 


The   value   of   w       and   uo        were    substituted    into   Eq .     (3.4)    and   yielded 
kl  k2 


0 

M,= 


I    *    ~—  (3-13) 


Where   ' 


M,     =  M        =  M      ,    recalling     CU&-    ~\i  fo  -  c£    „       Loci   of   constant   values 
of  M     were   plotted    in   relation    to   "a"    and   "b" .       (Fig.    3-8) 

K 

3. 7      N   Calculations 
From  Eq,     (3.  3) 


from   which, 

A/  =  -  a, dan  [    /**"  %  ]  ( 3. 14) 

Equation  (3.14)  was  then  the  basic  N  equation  for  the  pure  second  order 
system. 
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3.8   N  Determination 
P        

0)  was  previously  calculated  by  Eq .  (3.6).   The  substitution  of 

Eq  .  (3.6)  into  Eq  .  (3.14)  resulted  in  the  solution  of  N   related 

directly  to  "a"  and  "b",  the  closed  loop  poles  of  the  system. 

:  —  a  ret  an 


N, 


h  -a 


a 


(3.15) 


Loci  of  constant  values  of  N  were  then  determined  for  all  possible 

P 
"a"  and  "b"  values.   (Fig.  3-9) 

It  is  pertinent  to  note  that  the  equations  presented  were  in  most 
cases  not  solved  as  stated..   Algebraic  manipulations  were  performed  to 
make  the  equation  compatible  to  digital  calculations  by  writing  the 
equations  in  terms  of  "b"  as  a  function  of  "a"  and  the  particular 
parameter . 


3.9   N   Determination 
U)n 

N   was  determined  in  a  manner  similar  to  N 
am  p 

terms  of  "a"  and  "b"  was  stated  in  Eq .  (3.10).   This  relationship  of 
a  was  substituted  in  Eq .  (3.14)  and  resulted  in 


The  value  of  o>   in 
n 


ton 


da  n 


z  cii  &■+&*- 


a+b-  faar+b  J  J 


(3.16) 


which   simplified   to 


i  2. 


-arctan  ■ ,  2a.Va-i-.fc' 
L         0 


■ 


or 

N        =   -90° 
ion 

Thus,  N    for  a  pure  second  order  system  was  a  constant  and  equal  to  -90 

am       r  J 
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3.10   N,   Determination 

bw 

The  equation  for  N,   in  terms  of  "a"  and  "b"  was  obtained  by 

bw  J 

substituting  Eq .  (3.8)  into  Eq .  (3.14).   This  resulted  in 


(3.17) 


(Fig,   3-10) 


3.11   N  Determination 


Equalion  (3.12)  determined  the  values  of  CD  when  M  ■  1.0  for  non-zero 
uencies.   When  the  expressioi 
was  substituted  into  Eq .  (3.14), 


frequencies.   When  the  expression  of  uo  ,  determined  from  Eq .  (3.12), 


'l 


(3.18) 


Solutions  to  Eq .  (3.18)  were  plotted  as  constant  N1  loci  for  all  values 
of  "a"  and  "b".  (Fig.  3-11) 


3.12   N,  Determination 
k 

As  stated  above,  the  k  values  selected  were  .5  and  1.5  where 

£0,  — :"^|<  CO jj,  ,     In  a  manner  similar  to  the  determination  of  N,  above, 


(3.19) 


from  which  N  values  were  determined  explicitly  in  terms  of  "a"  and  "b" 
for  both  values  of  k.   (Figs.  3-12  and  3-13) 


3.13  Use  of  the  Curves 

At  this  point,  all  the  selected  parameters  necessary  to  sketch  the 
M  and  N  contours  have  been  determined  from  the  location  of  the  closed 
loop  poles.   With  the  values  of  the  closed  loop  poles  known  and  the 
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curves  plotted  as  described  above,  a  rapid  graphical  approximation  of 

the  frequency  response  of  the  pure  second  order  system  can  be  constructed. 

For  convenience,  the  following  recommended  procedure  is  presented  for 

the  use  of  the  pure  second  order  curves: 

a.   Enter  Fig.  3-2  with  the  complex  pole  coordinates  -a  +  jb,  and 

obtain  M    .      Enter   Fig.    3-3  with   -a  +   jb,    and  obtain   uu    .       (Note    that   OJ 

P  P  P 

is    the   ordinate    in   Fig.    3-3.) 

b        Similarly  obtain   ID.        (M  =    .707)    from  Fig.    3-4,    U)     and  M 

bw  n      oon 

from  Figs   3-5  and  3-6  respectively. 

c.  Obtain  uo   (M  =1.0)  from  Fig.  3-7,  entering  with  the  UJ 
determined  above, 

d.  Obtain  M,  for  both  k  =  .5  and  k„  =  1.5  from  Fig.  3-8,  entering 
with  "a"  and  "b".   Determine  oj   ,  (  u>^  -~^  ,5~  ^COp)9 

and        <v  (wk^'$^~'cop)> 

from  Fig.  3-7,  entering  with  the  value  of  UJ  obtained  from  step  (a)  above 

The  M  contour  for  the  specified  closed  loop  poles  may  then  be 

sketched,  having  easily  obtained  the  M  and  corresponding  OJ  for  six  known 

points  and  realizing  that  M  =  0  and  has  a  slope  of  zero  at  (l)  =  0,  and 

also  that  the  slope  is  zero  at  cjo  and  at  u)  =  °°. 

P 

The  N  contour  points  are  obtained  in  the  same  manner  by  using 

Figs.  3-9  through  3-12  to  obtain  the  values  of  N  corresponding  to  the 

same  frequencies  found  for  the  M  contour  above.   The  one  exception  is 

N   ,  which  is  always  equal  to  -90   for  this  system.   Using  these  N, 
lOn 

U)  values  and  knowing  that  N  =  0  and  has  zero  slope  for  uo  =  0 ,  and  also 
that  N  =  -180   for  U)  =  °°  for  a  pure  second  order  system,  the  N  contour 
can  now  be  accurately  sketched. 


Note  that  this  paper  does  not  present  curves  which  will  give  any 

values  at  frequencies  higher  than  (D   ,  since  it  is  unlikely  that  there 

bw 

would  ever  be  any  interest  in  the  region  beyond  the  bandwidth  frequency, 
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3  14  Example 

Given:    -a  +  jb  =  -1  +  j2 

Sol ution: 
SPECIFIC  DATA 


TABLE  OF  PARAMETER  VALUES 


X 

0) 
X 

(JU     Fig.    No. 

X 

M 

X 

M      Fig,    No. 

X 

N 

X 

N      Fig.    No. 

X 

p 

If! 

3-2 

1.Z5 

3-3 

-60° 

3-9 

bw 

Z$* 

3-4 

.707 

N.A. 

-J*3° 

3-10 

CD 

n 

z>zs 

3-5 

1.12 

3-6 

-90° 

N.A. 

1 

/■■: 

3-7 

1.0 

N.A. 

-AV 

3-11 

ki 

3-7 

/.//$- 

3-8 

-J^° 

3-12 

k2 

Z.IO 

3-7 

2.//5" 

3-8 

-a^° 

3-13 

GENERAL  DATA   (for  all  pure  second  order  systems) 


uu 

M 

dM/duo 

N 

dN/dou 

0 

CJU 

p 

00 

1.0 

M 
P 

0 

0 
0 

0 

0° 

0 

-180° 

0 

Fig   3-14  (page   17  )  compares  the  frequency  response  obtained  from 
this  rapid  approximation  method  to  the  actual  M  and  N  contours  obtained 
from  a  digital  computer  program  of  the  general  M  and  N  equations. 

Note :   Although  a  scale  is  shown  on  these  pure  second  order  curves  for 
convenience,  the  loci  are  really  independent  of  scale  as  they  depend 
only  on  the  ratio  of  "a"  to  "b"    Hence,  any  convenient  scale  for  O, 
jet),  and  CD  may  be  substituted  for  the  scales  shown  as  long  as  all 
scales  are  in  agreement 
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Pig.  3-6 
l^n   for  Pure  Second  Order  System 
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25 
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7^ 


Fig.  3-10 
N,   for  Pure  Second  Order  System 


m 
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Fig.  3-11 
N^  for  Pure  Second  Order  System 
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4.   THE  SECOND  ORDER  SYSTEM  WITH  ONE  ZERO 

The  addition  of  a  real  zero  to  a  pure  second  order  system 
resulted  in  a  system  whose  closed  loop  transfer  function  was 

~A — 7 \a~T  (  1) 


The  singularities  of  this  system  were  plotted  on  the  s-plane  in 

2    2 
Fig  4-1.   The  constant  a  +  b   ,  usually  interpreted  as  system 

d 
gain,  reflected  the  assumption  of  unity  feedback  as  stated  previously. 

By  letting  s  =  juu,  the  frequency  response  form  of  Eq.  (4.1) 
resulted  in 


t-cWrn — 7f — rr  (4.2) 


This  expanded  and  rearranged  became 


(4.3) 


r/<  \  -^  (diA <^ £) 4 j  (a? £-'  +  o>ba) 

3 
Dividing    the   numerator   and   denominator  of  Eq .     (4.3)    by   d    , 

F/.   ,     (d  VJw?+TdV  <44> 


Id1 


From  Eq.     (4,4),    it   was    seen   that  "a",  *b",  and  to  may  all   be   referred   to    the 
zero    location,    d. 
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s-plane 


■e- 


x- 


+  j"> 


-x^ 


</ 


k- 


x-- 


a 


Fig.  4-1    Poles  and  Zero  of  the  Second  Order  System  with  One  Zero 
(s-plane) 
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From  Eq.  (4.4) , 


™i*  TTT 

1   |d|  (4-5> 


•7/1-^ 


*"  \d\  (4.6) 


and  XI  — 


^  a) 


|d|  (4'7> 


These  quantities  were  illustrated  in  Fig.  4-2.   Equation  (4.4)  then 
became 


r-/.  \   On,  +7/ia)-h1tJQ./ft,  -f~TU^  J 

'ctpi-  ;    i — t     ^      ;         ^    ''  (4.8) 


Thus  the  frequency  response  for  the  second  order  system  with  one  zero 
became  a  function  of  the  "referred"  location  of  the  complex  roots  and 
the  "referred"  frequency.   This  means  then  that  the  frequency  response 
was  the  same  for  all  systems  having  the  same  relative  pole-zero  con- 
figuration when  the  frequency  was  also  referred  to  the  zero  location. 
(Fig.  4-3) 

4. 1   M  and  Q  Calculations 

Equation  (4.8)  was  expressed  in  terms  of  its  magnitude  and  phase 
angle  by 

Yc ( J4"  s  7~3 ''.■  w, "Y"  =  M €*       (4.9) 
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s-plane 


m„ 


a 

Ml 
h 


% 


k 


?>?i.d 


d 


rntd  — 


A    !<*> 


J 


-*-oJ 


m-plane 


r 


K mx 


X 


-1 


A  im 


Jm* 


*->fy 


Fig  4-2   Referred  Quantities  of  m  and  rop ,  Second  Order  System 
with  One  Zero 
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System  A 
a  =  1.0 
b  -  2,0 
d  =  1.0 


System  B: 
a  -  .01 

b  -  .02 
d  =  .01 


ml  =/d/=  1-° 


mi =  /ir  1-° 


m2  -  /?/=  2.0 


m, 


7a/-  2-° 


(JU 


0 

1 

2 

3 

4 

Q  (Q  -/l7) 

0 

1 

2 

3 

4 

u)     (d   -   1.0) 

0 

.01 

.02 

.03 

.04 

(DB    (d    =     .01) 

Fig.  4- 3  Frequency  Response  of  Two  Systems  with  the  Same  Relative 

Pole-Zero  Configuration,  Second  Order  System  with  One  Zero 
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From  this , 


A 


IVI  =  —7=        ■ ;       =r  (4.10) 


Equation  (4.10)  was  the  basic  M  equation  for  the  second  order  system 
with  one  zero. 


4.2   M  and  Q  Determination 
P      P 

The  equation  for  Q  was  obtained  by  first  squaring  Eq.  (4.10), 
then  differentiating  with  respect  to  &  and  setting  the  numerator  of 
this  equal  to  zero.   This  reduced  to 

jLf= Y-i  +Vw+i^)1-2.6nr-»*) +i '        (4-u) 


Substituting  this  into  Eq.  (4.10)  and  simplifying, 


"f>~~    \\         \     A  i     x  4  %  ==i  I  (4.12) 


Loci  of  constant  values  of  m„  were  solved  by  digital  computer  and  were 

plotted  against  m.  versus  M  and  Q  .   (Figs.  4-4  and  4-5). 
1         p      p 

4.3  Q,   Determination 
bw 


The  equation  for  the  referred  bandwith  frequency  was  obtained  by 
substituting  M  ■  .707  into  Eq.  (4  10)  and  solving  for  CI     .   After 
considerable  algebra,  this  reduced  to 

+  tmi^LJ.+  |]-4w>H  \     (4.n) 

(Fig.  4-6)  Z  J 
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4.4     MU     and  Q     Determination 
Oi n 

The    referred  undamped  natural    frequency  was   defined   for    this 

system  as 


Sln.=  ~i'h^7H 


■j 


(4.14) 


This  expressed  0  as  a  function  of  m  and  m_  .   (Fig.  4-7) 

Substituting  Eq .  (4.14)  into  Eq .  (4.10)  and  simplifying  yielded 


Kr-\lU+^+P"^'  (*.«) 


(Fig.  4-8) 


4.5   fi  Determination 

As  defined  previously,  Q  was  the  referred  frequency  other  than 
zero  at  which  M  is  unity.   The  equation  for  Ci  was  found  by  letting 
M  ■  1.0  in  Eq.  (4.10)  and  solving  for  0  . 


-ft  x  =  i  (rf+yntf-z(r>h-7yiZ)  (*■ 16) 

(Fig.  4-9) 

4.6  R  and  CI  Determination 

The  use  of  the  curves  mentioned  up  to  this  point  for  this  class 

of  system   provided  the  four  main  points  required  to  sketch  roughly  the 

M  contour.   In  order  to  "fill-in"  the  blank,  portions  of  the  contour  on 

both  sides  of  d   ,  it  was  decided  to  provide  curves  for  three  additional 

P 
points  on  the  M  contour  in  the  same  manner  as  was  done  for  the  pure 

second  order  system.   The  use  of  these  additional  points  was  not 

necessary  to  sketch  the  frequency  response;  however,  their  use  provided 

a  mere  accurate  approximation  of  the  M  contour 
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As  in  the  pure  second  order  system,  the  following  definition  was  made 

where  k  is  a    constant 

Substituting  into  Eq .  (4  11)  for  Q  , 


(4.17) 


IU  =  -V-I<  +  ky(7n?+7niy~z(7n?-mZ)+  1 


(4.18) 


Values  of  .5,  1.5  and  2.5  were  selected  for  k.   M,  solutions  were 

k 

obtained  for  each  k  value  by  solving  Eq .  (4.18)  for  £i  and  then  sub- 
stituting this  value  into  the  basic  M  equation  with  the  same  values  of 
m  and  m,-   This  was  accomplished  handily  by  one  digital  program. 
(Figs   4-10,  4-11,  4-12,  and  4-13) 

4, 7   N  Calculations 

Referring  to  the  basic  closed  loop  transfer  function  in  frequency 
form,  Eq.  (4.9),  the  general  phase  angle  equation  for  the  second  order 
system  with  one  zero,  was  written  as 


N  ~  '    " 


(4.19) 


The  previously  obtained  parameter  frequency  equations  (ft  ,  ft,  ,  ft  , 

p    bw   n 

ft  ,  and  the  three  ft,  's)  were  solved  by  digital  computer  programs.   The 
resulting  parameter  frequencies  were  then  substituted  into  the  general 


N  equation  (4,19)  along  with  the  corresponding  values  of  m  and  m  .   The 
following  figures  resulted,  all  plotted  with  constant  m.,  loci  against  m 
versus  Nr 


Fig,  4-14 

Fig.  4-15 

Fig.  4-16 

Fig,  4-17 

Fig  4-18 

Fig.  4-19 

Fig-  4-20 


N      as    function   of  mn    and  m« 
p  12 

N.       as    function   of  m,    and  m_ 
bw  12 


NU     as    function    of  m,    and  itk 
ifo  12 


N.    as 


function   of  m     and  m0 
1  Z 


N.       as    function   of  m:    and  m~  for   kn    =    .  5 
k^                                          12  1 

N,       as    function   of  mn    and  m0  for   k_    =    1.5 
k2                                          12  2 

N        as    function   of  m,    and  m.,  for   k  _.    =   2.5 
ko                                        12  3 
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4.8  Example 

Given:   Complex  poles;   a  +  jb  =  -2.0  +  j  2.0 

Real  zero;   -d  =-2.0 

Solution: 
SPECIFIC  DATA 

TABLE  OF  PARAMETER  VALUES 


X 

Q 

X 

0) 
X 

Q     Fig.    No. 

X 

M 

X 

M      Fig.    No. 

X 

N 

X 

N      Fig.    No. 

X 

p 

bw 

fi 
n 

2.9 
1.42 

2.2, 

Z84 

4-4 
4-6 
4-7 

123 
.707 

123 

4-5 

N.A. 

4-8 

-23  J  ° 
-6/.00 
-3S.2* 

4-14 
4-15 
4-16 

1 

2.0 

is 6 

4-9 

4-10 

1.0 

1.22 

N.A. 
4-11 

-20.5° 

4-17 
4-18 

k2 

1.2S 

Z/o 

4-10 

1-24 

4-12 

-  3F.2° 

4-19 

k3 

i>7* 

3-t8 

4-10 

t.D?S 

4-13 

-*7.e° 

4-20 

Obtain  O)  by  (JO  =  /d/  C 
x     x        : 


GENERAL  DATA   (for  all  Second  Order  plus  One  Zero  Systems) 


Fig   4-21  (page  39  )  compares  the  frequency  response  obtained  from 
this  rapid  approximation  method  to  the  actual  M  and  N  contours 
obtained  from  a  digital  computer  solution  of  the  general  M  and  N 
equations . 
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Mp  for  Second  Order  System  with  One  Zero 
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Fig.  4-8 
M^r  for  Second  Order  System  with  One  Zero 
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Pig.  I4.-9.  -fL1  for  Second  Order  System  with  One  Zero 
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Fig,  lj.-12.   Mk2  f°r  Second  Order  System  with  One  Zero 
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Fig.  4-13.   Mi^o  for  Second  Order  System  with  One  Zero 
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Nbw  for  Second  Order  System  with  One  Zero 
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f_^n  for  Second  Order  System  with  One  Zero 
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5        THE,    PURE   THIRD   ORDER   SYSTEM 

The  addition  of  a  real  root  to  the  pure  second  order  system 
resulted  in  a  new  system  which  will  be  referred  to  as  the  pure  third 
order  system.   The  closed  loop  transfer  function  of  the  pure  third 
order  system  had  one  pair  of  complex  conjugate  poles  and  one  real  pole, 

The  closed  loop  transfer  function  for  the  unity  feedback  third 
order  system  was  written  as 


• 


(5.1) 


Fig   5-1  defined  these  roots  on  the  s-plane. 
Equation  (5.1)  was  simplified  as  follows: 


/  c  I  >     ~~    — — '  (  5  2^ 

.  -«•  •    -tb  J 

and  expanded  became 

b 


54      4  j  +  ?(cl  +  Z.;-  ■  4   )  +  c  "  .^N-i*)       (5-  3) 


Letting  s  =  jcu  in  Eq.  (5.3)  resulted  in  the  frequency  response  form 
of  the  transfer  function 


*     -  y  \  hi  !  (5-4> 


ISO 


or 


^ 


(a  "4 1 


'       '  J ''  •  ■  L°5}         (5'5) 
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Fig   5-1   Poles  of  the  Pure  Third  Order  System  (s-plane) 
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Numerator  and  denominator  of  Eq .  (5.5)  divided  by  c" 


r  (v  o 


-a 

cU'      c     cv    c3 


(5.6) 


Observing  that  a  power  of  c  exists  in  all  terms  corresponding  to 
the  same  power  of  "a",  "b"  and  0);  the  following  definitions  were 
made: 


71 


a  b 


(5.7) 


(5.8) 


M 

O  A  a) 

Icj  (5-9) 

These  quantities  were  illustrated  in  Fig.  5-2.   Equation  (5.6)  then 
became 


I  J -J)- 


[n*-H(l~Sl(  n,+]    -kjUl (?  i  Vm  ,  4  :>  ,""  S  -J}5 


^  M^! 


IN 


(5.10) 


n1  and  n„  were  the  complex  pole  locations  referred  to  the  real  pole 
location  and  ft  was  the  referred  frequency.   The  frequency  response 
for  the  pure  third  order  system  was  then  a  function  of  the  referred 
complex  pole  and  the  referred  frequency.   Thus  all  systems  having  the 
same  relative  pole  configuration  will  have  identical  response  if  the 
frequency  were  also  referred  to  the  real  pole  location. 


5. 1  M  and  Q   Calculations 

The  magnitude  portion  of  Eq .  (5.10),  after  algebraic  manipulation, 
reduced  to 


M  = 


3- 


7li  +7li 


V*       --..Mi 


(5.11) 
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Fig.  5-2  Referred  Quantities  n..  and  n2 ,  Pure  Third  Order  System 
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Equation  (5.11)  was  the  basic  M  equation  for  the  pure  third  order 
system. 

5.2   Discussion  of  the  M  Contours 

It  was  found  that  the  M  contour  for  a  pure  third  order  system 

will  be  one  of  three  basic  "shapes"  as  shown  in  Fig.  5-3a,  b  and  c, 

depending  on  the  values  of  n1  and  n„ .    Figure  5- 3a  was  the  situation 

when  there  was  no  peak.   No  value  of  Q     was  obtained  and  hence  no  values 

P 
for  Cl    or  the  Cl  '  s.   Figure  5-3b  was  similar  to  the  second  order  M 

contour  shape  and  had  a  single  peak.   For  this  general  "shape"  points 

were  obtained  for  CJ  ,  Cl  ,  Q   ,  Q   ,  Cl  tiL     ,  and  Cl  . 

The  third  general  "shape"  is  Fig.  5-3c,  which  resulted  from  the 

Q  equation  (5.12)  which  yielded  two  non-zero  frequencies  at  which 

the  M  contour  slope  was  zero,  one  being  a  "peak"  frequency  and  the 

other  being  a  "minimum"  frequency.   This  latter  frequency  was  designated 

Q     to  avoid  being  confused  with  the  Cl   .   It  was  also  discovered  that 
m  p 

with  this  shape  of  M  contour,  it  was  possible  to  have  M  less  than 
unity.   Another  possibility  was  to  have  M  less  than  M,   and  at  the  same 
time  M  greater  than  one.   This  resulted  in  two  values  of  Q  ,  two  values 
of  Cl  ,  M,  ,  Cl  ,  and  M,   (one  relative  to  M  and  the  other  relative  to 


u      ,  M,  ,  Cl  ,  and  M,   (one  relative  to  M 


M_) ,  and  three  values  of  Cl   (one  for  each  crossing  of  the  M  =  .707  line), 


m 


in  addition  to  the  usual  points  at  Q   ,  C   ,  and  Q  . 

p   m       n 

Figure  5-4  might  be  thought  of  as  the  complex  n-plane  and  shows  the 

regions  in  which  certain  n..  and  n_  combinations  produce  one  of  the 

three  above  M  contours.   This  figure  indicates  the  general  shape  of  the 

M  contours  that  result  from  the  known  n1  and  n„  values. 

5.3  M  and  Q  Determination 
_J2 2 

The  equation  for  Ci  was  obtained  by  squaring  Eq .  (5.11),  then 
differentiating  with  respect  to  Q  and  setting  the  numerator  of  this 
equal  to  zero.   This  reduced  to 
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Fig.  5-3  The  Three  Basic  M  Curve  Shapes  for  Pure  Third  Order  System 


(a) 


(b) 


(c) 


/*l      J 
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The  plus  and  minus  values  of  the  inner  square  root  quantity  were  both 

valid  and  resulted  in  two  possible  solutions  for  Ci   .   The  consequences 

P 
of  two  solutions  were  discussed  in  detail  in  paragraph  5.2  above. 

Loci  of  constant  values  of  n~  were  solved  by  a  digital  computer 

program  and  were  plotted  against  n.  versus  Ci     and  Ci     in  Figs.  5-5a 

1         p      m 

and  5- 5b    In  the  same  computer  program,  values  of  M  were  obtained 

for  corresponding  values  of  n  ,  Ci   ,   and  n_  by  substituting  these 

known  quantities  into  the  basic  M  equation,  Eq .  (5-11).   Figures  5-6a, 

5-6b.  and  5-6c  contain  M  plotted  as  a  function  of  n,  and  n~ .   No 

P  12 

specific  M  equation  was  derived  due  to  the  obvious  complexity 
involved  in  substituting  Eq .  (5.12)  into  Eq .  (5.11)  for  Ci. 


5.4  CI       Determination 

bw 

The  equation  for  the  referred  bandwidth  frequency  was  obtained 
by  substituting  M  =  .707  into  Eq   (5„11)  and  solving  for  n  as  a  fur 
of  Cl       and  n_   This  involved  extensive  algebra  which  reduced  to 


Ti,= 


%\  4\iA  ")-%L 


•%»■  *  A^e2.^4l)-f-Q-hc.,-  (*3.-2*  J  ~  ft-oH   L        (5.13) 


w' 


The  +  signs  were  both  valid  here  also  and  their  effect  will  be  discussed 
later.   (Fig.  5-7) 

5.5  M^  and  Cl     Determination 
L2n      n 

The  undamped  natural  frequency  for  the  pure  third  order  system  was 
defined  as 


(5.14) 


Cl 


!  was  expressed  as  a  function  of  n.  and  n,.   (Fig.  5-8) 
n  \  I 


6.3- 


Substituting  this  into  Eq .  (5.11)  and  simplifying  yielded 


(5.15) 


V  4**6*N7£  )(n?+nl  •+  l) 


Mr  was  thus  expressed  as  a  function  of  n.  and  n0 .   (Figs.  5-9a  and  5-9b) 
i<n  1      L 

5.6  Cl     Determination 

The  referred  frequency,  other  than  zero,  at  which  M  ■  1.0  was 
defined  as  Q.  . 

Letting  M  ■  1.0  in  Eq.  (5.11)  and  solving  for  fi.  resulted  in 


J) 


I  * y-nX-.5  ±inl~rf-4rfnl+wZ5  (5. 16) 


(Figs.  5-10  a,  b,  c,  and  d) 
5.7  M,  and  CI  Determination 


For  the  same  reasons  as  discussed  in  Section  2.4,  it  was  found 
rable  to  obtain  ' 
was  done  by  defining 


desirable  to  obtain  "fill-in"  points  on  both  sides  of  Cl   .   This 

P 


ft  -  kQp  (5.17) 

which  eliminated  the  square  root  of  k  as  used  for  the  previous  two 
classes  of  systems. 

The  values  of  k  chosen  were  .8,  1.4,  and  1.8;  therefore 

CL      ■  8  Qp  (5.18) 

kl 

Cl   ■  1-4  Op  (5.19) 

k2 

fi   -  1.8  C*p  (5.20) 

k3 
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Fig.  5-11  was  the  almost  trivial  plot  of  these  three  equations, 
included  for  optional  use  in  determining  the  £1  values  once  Q   has 
been  determined. 

Loci  for  the  three  M,  's  were  calculated  by  a  digital  program  which 
utilized  Eqs.  (5.11),  (5.12),  (5.18),  (5.19),  and  (5.20);  thus  it  was 
not  necessary  to  attempt  to  obtain  any  specific  M,  equations. 
(Figs   5-12  a,b,  5-13  a,b,  5-14) 

5.8  Determination  of  N  Curves 

Referring  to  the  basic  closed  loop  transfer  function  in  frequency 
form,  Eq .  (5.10),  the  general  phase  angle  equation  for  the  pure  third 
order  system  was  written  as 


(5.21) 


The  previously  obtained  parameter  frequency  equations  (Q   ,  Ci     , 
Cl   ,    and  Q  )  were  solved  by  digital  computer  programs.   The  resulting 

X.  K 

parameter  frequencies  were  then  substituted  into  the  general  N  equation, 
Eq   (5.21),  along  with  the  corresponding  values  of  n.  and  n~  .   The 
following  figures  resulted,  all  plotted  with  constant  n?  loci  against 
n1  versus  N: 

Fig-  5-15  a    N  as  a  function  of  n,  and  n0 

°  p  12 

Fig.    5-15  b  N     as   a    function  of  n,    and  n_ 

m  1  L 

Fig.     5-16  a  N.       as   a    function  of  n.    and  n„ 

bw  12 

Fig.    5   16  b  N,       (extended)   as   a    function  of  n.    and  n0 

bw  12 

Fig      5-17  a  N.    as  a    function  of  n1    and  n. 

Fig.    5-17   b  N      (expanded)    as   a    function   of   n      and  n_ 

Fig.    5-17   c  N       (2nd  value) as   a    function  of  n1    and  n„ 

Fig      5-18  a  Nb      for  maximum    (ft  ) 

k1  P 
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Fig.    5-18   b  N.       for  minimum    (Q  ) 

K,  m 

Fie.    5-19   a  Nl      for  maximum    (Cl  ) 

k2  p/ 

Fig.    5-19    b  N,       for  minimum    (Q  ) 

6  k  nr 

Fig.    5-20  N,       for  maximum    (Q  ) 

^3  P 

An  equation  for  N^s  was  obtained  by  substituting  Eq.  (5.14)  into 
the  basic  N  equation,  Eq .  (5.21), 

NJln  =- arc-ban [~^J 

—  —  o^rcian  I  — — : —I  (5.22) 

Loci  of  Nq^  were  plotted  as  a  function  of  n.  and  n~  .   (Fig.  5-21) 
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5. 9   Example 

Given:   Complex  poles- 


Real  pole: 


Solution : 
SPECIFIC  DATA 


■a  +  jb  =  -.8  +  j  4.0 
-c   =-  2.0 


nl  "  Id 


n2  "  Id 


=  .4 


-  2.0 


TABLE  OF  PARAMETER  VALUES 


X 

X 

(JO      k  irk 

X 

Q     Fig. No. 

X 

M 

X 

M      Fig. No. 

X 

N 

X 

N      Fig. No. 

X 

p 

/  86 

s./6 

5- 5a 
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5-6a,b 
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5- 15a 
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A  #6 
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5-15b 

bw 
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bw 

N   A 
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.707 

N.A. 

N.A. 

— ■ — 

bw 

N.A. 



.  707 

N.A. 

N.A. 

n 

Z.a* 

4.o & 

5-8 

/■  t2 

5-9b 

-JS"4 

5-21 

1 

4>2i 

5-10a,b 

1.0 

N.A. 

-166  :' 

5-17a,b 

1 

Mf 

Z.JO 
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1.0 

N.A. 

-S 

5-17c 
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5-11 
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5-12a 

-8?* 

5-18a 
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0.76 
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5-11 

9/2 

5-12b 
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5-11 
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GENERAL  DATA  (for  all  pure  third  order) 


u)   (or  Q) 

M 

dM/duo 

N 

dN/doo 

0 

1.0 

0 

0° 
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M 
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00 
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-270° 

0 

Fig.  5-22  (page   69  )  compares  the  frequency  response  obtained  from 
this  rapid  approximation  method  to  the  actual  M  and  N  contours 
obtained  from  a  digital  computer  program  of  the  general  M  and  N 
equations . 
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Fig.  5-22 
Example (Frequency  Response 
for  Pure  Third  Order  System ) 
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Fig.  5-6a 
M^  for  Pure  Third  Order  Syst 
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Fig.  5-6c  • 
Mja  for  Pure  Third  Order  System 
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-       Fig-  5-7 
JLbw  for  Pure  Third  Order  System 
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Fig.  5-9a 
M_Q.n  for  Pure  Third  Order  System 
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Fig.  5-10a.  SL^   for  Pure  Third  Order  System 
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Pig.    5-10c.   jfL-^   for  Pure   Third   Order   System 
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Fig.    5-12a 
Mu    for  Pure   Third  Order  System 
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Fig.  5-l6b 
for  Pure  Third  Order  System 
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Nj_  for  Pur©  Third  Order  System 


Fig.   5-17 b 
Pure  Third  Order 
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Fig,  5-19a 
Nk2  for  Pure  Third  Order  System 
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Fig.  £-19b 
Nk5  for  Pure  Third  Order  System 
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Fig.   5-21 
NJln  for  P^©   Third  Order  System 
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6,   CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  conclusions  were  made  based  upon  the  investigation 
of  the  three  classes  of  servomechanisms  studied  in  this  thesis: 

a.  All  pure  second  order  servomechanisms  with  the  same  relative 
complex  pole  location  will  have  the  same  frequency  response. 

b.  For  the  pure  second  order  servomechanism,  the  parameter  loci 
plotted  are  made  applicable  to  any  location  of  complex  poles  in  the 
left  half  of  the  s~plane  by  changing  the  scales  on  the  complex  plane 
and  on  the  parameter  frequency  loci. 

c.  The  frequency  response  for  servomechanisms  with  three 
singularities  can  be  expressed  as  a  function  of  the  referred  real  and 
imaginary  parts  of  the  complex  poles  and  the  referred  frequency  where: 

1   The  referred  real  part  of  the  complex  poles  is  the 
actual  real  part  of  the  complex  poles  divided  by  the 
magnitude  of  the  real  pole  or  zero. 

2.  The  referred  imaginary  part  of  the  complex  poles  is  the 
positive  imaginary  part  of  the  actual  complex  pole 
divided  by  the  magnitude  of  the  real  pole  or  zero. 

3.  The  referred  frequency  is  the  actual  frequency  divided 
by  the  magnitude  of  the  real  pole  or  zero. 

d.  All  servomechanisms  with  three  singularities  which  have  the 
same  relative  pole -zero  locations  will  have  the  same  frequency  response 
relative  to  the  referred  frequency. 

e.  Curves  for  determining  values  of  several  optional  "fill-in" 

parameters  were  included  in  addition  to  the  normally  accepted  parameters 

of  M   ,  U)   and  U),  . 
pur  p      bw 

f   An  accurate  graphical  approximation  of  the  frequency  response 

of  a  servoniechanism  can  be  rapidly  obtained  using  the  curves  produced 

in  this  thesis, 

It  is  recommended  that  the  following  areas  be  considered  for  further 
investigation  using  the  technique  presented  in  this  paper; 

a.   Servomechanisms  with  more  than  three  singularities,  and 
b..   Systems  other  than  servomechanisms. 
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8 .   Appendix  A 

The  Control  Data  Corporation  Model  1604  digital  computer 
important  characteristics  are  summarized  as  follows: 

1.  All  transistorized 

2,  48  bit  word  length 
3o   Stored  program 

4    Single-address  logic,  two  instructions  per  48-bit  word 

5.   Indirect  addressing  feature 

6    32,768   48-bit  words  of  magnetic  core  storage: 

(a)  Stored  in  two  independent  16,384  word  banks, 
alternate  phased 

(b)  4.8  micro  -  seconds  effective  cycle  time 
(representative  program) . 

(c)  6,4  microseconds  total  cycle  time. 
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